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ABSTRACT 



Aims. We studied the physical properties of three clusters of galaxies, selected from a BeppoSAX wide field camera (WFC) survey. 
These sources are identified as 1RXS J153934.7-833535, 1RXS J160147.6-754507, and 1RXS J081232.3-571423 in the ROSAT 
All-Sky Survey catalogue. We obtained XMM-Newton follow-up observations for these three clusters. 

Methods. We fitted single and multi-temperature models to spectra obtained with the EPIC-pn camera to determine the temperature, 
the chemical composition of the gas and their radial distribution. Since two observations were contaminated by a high soft-proton 
background, we developed a new method to estimate the effect of this background on the data. 

Results. We present the temperature and iron abundance of two of these three clusters for the first time. The iron abundance of 
1RXS J153934.7-833535 decreases with radius. The fits to the XMM-Newton and Chandra data show that the radial temperature 
profile within 3' towards the centre either flattens or lowers. A Chandra image of the source suggests the presence of X-ray cavities. 
The gas properties in 1RXS J160147. 6-754507 are consistent with a flat radial distribution of iron and temperature within 2' from 
the centre. 1RXS J081232. 3-571423 is a relatively cool cluster with a temperature of about 3 keV. 

Conclusions. The radial temperature and iron profiles suggest that 1RXS J153934.7-833535 is a cool core cluster. The Chandra 
image shows substructure which points toward active galactic nuclei (AGN) feedback in the core. The flat radial profiles of the 
temperature and iron abundance in 1RXS J160147. 6-754507 are similar to the profiles of non-cool-core clusters. 

Key words. X-rays: galaxies: clusters - Galaxies: clusters: individual: 1RXS J081232.3-571423,1RXS J153934.7-833535, 
1RXS J160147. 6-754507 - Galaxies: abundances - Galaxies: clusters: intracluster medium 



1. Introduction 

Clusters of galaxies are the largest gravitationally bound ob- 
jects in the universe. They are the result of the growth of over- 
densities that were present in the early universe. About 80% of 
the mass in clusters consists of dark matter. Until now it has been 
impossible to directly measure it. The baryonic matter in clusters 
of galaxies resides in the stars and a hot diffuse plasma (hereafter 
referred to as gas). The stars account for only 1-3% of the total 
cluster mass, while the gas is responsible for 10-15% (e.g. David 
et al. 1990; Allen et al. 2002). This gas, called the intra-cluster 
medium (ICM), has temperatures of ~ 10 7_8 K, high enough to 
emit X-ray radiation. 

Originally, clusters of galaxies were discovered based on 
their optical properties (e.g. Abell 1958). In the 1960s and 1970s, 
X-ray emission was detected from positions or directions coin- 
ciding with those of known clusters of galaxies (Byram et al. 
1966; Bradt et al. 1967; Kellogg et al. 1971). Cavaliere et al. 
(1971) stated that this emission was indeed originating from the 
clusters. In the past ten years, the instruments on board of the 
XMM-Newton and the Chandra X-ray telescopes enabled the 
study of the physical and chemical properties of the gas in more 
detail. 

In the X-ray band, two general catalogues have been com- 
piled. In the hard X-ray band (0.5 to 25 keV), HEAO-1 surveyed 
the sky with a detection threshold of 5 x 10~ 15 WtrT 2 (Wood 
et al. 1984). The soft X-ray band (0.1 - 2.4 keV) was probed in 
the ROSAT All Sky Survey, with a flux limit of 5 x 10~ 16 W trT 2 
(Voges et al. 1999). The ROSAT observations yielded various 



catalogues of clusters of galaxies (e.g. Reiprich & Bohringer 
2002; Ebeling et al. 2002; Bohringer et al. 2004). 

We present the results of a spectral analysis of three clus- 
ters of galaxies, selected from observations with the wide field 
cameras (WFCs) on board BeppoSAX. The Italian-Dutch X-ray 
observatory BeppoSAX (Boella et al. 1997) operated between 
June 1996 and April 2002 and carried two WFCs (Jager et al. 
1997). They had a field of view of 40° x 40° (full width at zero 
response), angular resolution of 5' (FWHM) and pointed in op- 
posite directions. They were sensitive in a bandpass of 2 to 28 
keV with a spectral resolution of 1.2keV at 6 keV (FWHM). 
The WFCs observed the whole sky. The exposure coverage was 
somewhat inhomogeneous. However, most positions on the sky 
have been observed for at least 1 million seconds. The sensitiv- 
ity level reaches 3 x 10~ 15 WtrT 2 . Therefore only the brightest 
sources of ROSAT are detected by the WFC. 

All data of this instrument were combined on time scales of 
six months as well as the whole mission and in energy ranges 
of 2-5, 5-9 and 9-28 keV. The combined data yielded detec- 
tions of a few of cataclysmic variables and tens of active galac- 
tic nuclei (AGN) and clusters of galaxies (see e.g. Verrecchia 
et al. 2007, for a catalogue based on part of this survey). Six 
objects were found to be ill-understood despite their brightness. 
The 2-10keV fluxes of these six sources lie in the 7 x 10~ 15 to 
3 x 10- 14 WirT 2 range. 

We followed up these bright sources using XMM-Newton. 
For one source we analysed the available Chandra archival data. 
Three of these six are clusters of galaxies. The physical prop- 
erties of the gas in these clusters were found using the X-ray 
spectra. 
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Table 1. Previously known parameters for the studied sources. 






xpix 

Fig. 1. Contours of an archival ROSAT HRI observation of J 1601. The 
three dots plotted are optically bright foreground stars. The observation 
is binned using a Gaussian with a 6" standard deviation. The nth contour 
corresponds to 0.8" times the maximal count rate. 



Throughout this paper, we use Hq = 70 km s _1 Mpc~', Qm = 
0.3, and Q.^ = 0.7. The abundances are in units of the proto- 
solar abundances by Lodders (2003). Unless specified otherwise, 
quoted statistical uncertainties are at a lcr (68%) confidence 
level. 



2. The clusters discussed in this paper 

Here, we introduce the three observed clusters and give an 
overview of the catalogues in which they appear. We refer to 
the three clusters by the first four numerals from their 1RXS- 
catalogue name preceded by the letter "J". The fluxes and lumi- 
nosities found by previous authors are quoted in Table 1 . 



2.1. J 1539 (1RXSJ1 53934.7-833535) 

This source is part of the REFLEX catalogue 
(RXCJ1539.5-8335, Bohringer et al. 2004), a cluster 
survey based on ROSAT data. 



2.2. J1601 (1RXSJ1 60147.6-754507) 

J 1601 is part of the ROSAT catalogue of clusters in the "zone 
of avoidance", meaning a location within 20° of the Galactic 
plane (CTZAJ1601. 7-7544, Ebeling et al. 2002). That is also 
the first detection and identification of this source as a cluster 
of galaxies. The source also appears in the RXTE All-Sky Slew 
Survey catalogue, (XSS J16019-7548, Revnivtsev et al. 2004). 
Figure 1 shows intensity contours of an archival observation of 
this source, using the HRI aboard ROSAT. 



2.3. J0812 (1RXSJ081232.3-571423) 

J0812 was first detected by the Einstein X-ray satellite as 
part of the Einstein Galactic plane survey (IE 081 1.5-5704, 
Hertz & Grindlay 1984). Furthermore, it appears in the Einstein 



Object 


WFC Flux " 


ROSAT flux" 


redshift c 


J1539 


0.82 


1.68 "•<• 


0.0728 d 






2.49 / 




J1601 


1.28 


2.41 f 


0.1530 * 


J0812 


1.14 


3.93 / 


0.0620 * 



" 2-10 keV, 10 14 Wm 2 for Crab-like spectrum (i.e. a power law with 

a photon index of 2) 

h 0. 1-2.4 keV, 10~ 14 WirT 2 

c Determined using optical spectroscopy 

d From REFLEX (Bohringer et al. 2004) 

e Using a Raymond-Smith model (Raymond & Smith 1977) 

/ From the RASS-BSC (Voges et al. 1999), assuming a power-law 

spectrum. 

* From CIZA (Ebeling et al. 2002) 



Table 2. Exposures of the XMM follow-up of the WFC data. 



Object 


time (ks) " 


filtered time (ks) 


date (Y-M-D) 


J1539 


17.7 


13.8 


2008-04-01 


J1601* 


9.12 


6.96 


2008-02-21 


J0812 


1.53 


1.53 


2008-03-01 



" Sum of good time intervals for central CCD. 

6 In our observation, the pointing was towards another unidentified 
ROSAT source (1RXS 1160107.6-754534) within the uncertainty 
circle of the WFC. 



Slew Survey (1ES 081 1.4-5704, Elvis et al. 1992) and is 
also part of the 2E catalogue (2E 081 1.4-5704, Harris et al. 
1990). The source was identified as a cluster of galaxies by 
Ebeling et al. (2002) and incorporated in the CIZA catalogue 
(CIZAJ0812.5-5714). 

3. Data reduction and analysis 

3.1. XMM-Newton EPIC 

A log of the XMM-Newton European photon imaging camera 
(EPIC; Striider et al. 2001) observations is given in Table 2. For 
all our observations the "medium" filter was put in place. For 
each observed source, we extract a lightcurve to look for soft- 
proton flares. Since emission from clusters of galaxies is not time 
dependent, we use the full energy range of the detector to extract 
this light curve. We compared the lightcurves within 3' around 
the cluster centre to the lightcurve of the remainder of the field. 
Both light curves show flaring at the same time scales which 
excludes the possibility of a variable point source affecting the 
light curve. For two observations (J 1539 and Jl 601), the residual 
soft-proton background was very high (see section 3.1.2). We 
excluded data during the flares on top of this high background. 
The original and filtered observation times are listed in Table 2. 

3.1 .1 . Standard data treatment 

The data reduction is done with version 8.0.1 of the XMM- 
Newton science analysis software (SAS). We used the calibration 
version of April 1st of 2009 and the standard EPIC-pn pipeline 
to create event files. During the observations, EPIC-MOS was in 
small window mode. We did not use these data because in this 
mode the source is larger than the field of view. We removed the 
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point sources by eye using a circle with a 15" radius around 
them. The next step involved correcting for the vignetting of 
the photons and protons. For the observation of J0812, which 
is not dominated by a high soft-proton background, we used 
the evigweight task to take this effect into account. This task 
assigns a theoretical vignetting factor to each photon based on 
its energy and position on the CCD. For the other two obser- 
vations, we took care of the soft-proton vignetting as described 
in Sect. 3.1.2. We selected only the single and double events 
(PATTERN<4) that were inside the field of view (FLAG==Q). For 
each separate spectrum, we generated an ancillary response file 
(ARF) and a response matrix file (RMF) using a detector map. 
For J0812, we took the detector map to be flat (i.e. without any 
spatial distribution) because spatial effects are already taken care 
of by the vignetting correction of evigweight. For J1539 and 
J 1601, we made an image of the event file and used that as a 
detector map. 

For the spectral reduction, we rebinned the spectra, sampling 
the instrument FWHM by a factor of three, requiring at least 
20 counts per bin to ensure the applicability of x 2 statistics if 
needed, and checked that spectral features were not undersam- 
pled. 

The next step was the background treatment. Since the 
sources are relatively compact, we took the background from an 
annulus ranging from 9' to 12', centred on the telescope point- 
ing. This is the outer edge of the field of view where we did not 
expect any residual cluster emission. We scaled the background 
spectrum with the ratio of the number of illuminated pixels in 
the observations to the background annulus. We subtracted the 
scaled background spectrum from the observation. Since there is 
a strong instrumental copper line at 8.04 keV in the background 
region, we ignored the energy range from 7.8 keV to 8.4 keV for 
analysis of the spectra to ensure that no radiation of this line was 
still present. For other instrumental lines, this effect is expected 
to be negligible in the energy range of the fit due to the statistics 
of the observations. 



3.1 .2. Special treatment of flared events 

Normally the intervals of soft-proton flaring in the light curve 
are excluded from the analysis. The observations of J1539 and 
J 1601 have a very high background. Inspection of the light 
curves of these observations for the region outside the 3' around 
the emission peak of the observation yields an average count rate 
of ~6c s -1 . In the SAS User Guide, the rule of thumb is to ig- 
nore all data for which the count rate is higher than 1 c s _1 . This 
would result in discarding all the data despite a positive detection 
of the source. In order to be able to use the data, we developed a 
method to estimate this background contribution. 

Kuntz & Snowden (2008) studied the soft protons in XMM- 
Newton data. They found a vignetting effect for soft protons in 
the EPIC-MOS detector. For the EPIC-pn, such a study has not 
yet been performed. Since the vignetting is caused by focussing 
of the soft-proton events by the mirror, we expect a vignetting 
similar to EPIC-MOS. Protons are vignetted less than X-ray 
photons. Therefore, standard SAS tasks for correcting the X-ray 
photon vignetting are not applicable. To be able to fit the spectra, 
we derived a vignetting factor for the soft-proton background of 
these data sets. We assumed that the X-ray photon background 
is negligible compared to the soft-proton background. Since the 
average background count rate is of a factor six higher than the 
maximum expected X-ray quiescent level, this is a reasonable 
assumption. Other simplifications are that the vignetting of soft- 
protons is taken to be energy-independent and that the vignetting 




0.8 0.9 1 1.1 

Background scale multiplication factor 



Fig. 2. Example of the Cash-statistic minimisation method, as used for 
cluster J1539. The multiplication factors are based on the region 3' 
around the cluster centre. 




Fig. 3. Chandra image of J1539, based on a 8ks observation using the 
ACIS instrument. There are no X-ray point sources in the region shown. 
The image is smoothed using a Gaussian with a standard deviation of 
2.5 pixels. The white circle represents the position of the central cD 
galaxy in this cluster. In the left panel, a hint of sub-structure is visible. 
The regions we used to fit the different parts of the spectrum are shown 
in the right panel. The projected radii of the elliptical "res." region are 
97kpc (width) and 89kpc (height). 



in the innermost 3' around the cluster centre is negligible. The 
radial behaviour of the soft-proton vignetting measured for MOS 
by Kuntz & Snowden (2008) is consistent with being flat within 
the inner 3'. 

In order to find the vignetting factor, we first scaled the back- 
ground to the region around 3' from the cluster centre, only 
based on the number of illuminated pixels. Then we created a set 
of spectra where we subtracted this background multiplied by a 
range of factors and fitted a model (as described in Sect. 3.3) to 
all of them. The C-statistic values of the fits as a function of the 
scale factor for J 1539 are shown in Fig. 2. We fitted a parabola to 
the data points to find the vignetting factor and the uncertainty 
on this value. The multiplication factor at the minimum of the 
parabola was taken as the best-fit value. The lcr uncertainty is 
then given by the scale factor value for which the C-statistic is 1 
above the minimal value (Cash 1979). 

3.2. Chandra 

Chandra data are available only for J1539. The 8ks ACIS ob- 
servation was taken in pointing mode on 24 June 2007. The 
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Chandra observation does not suffer from high background. This 
together with a high spatial resolution allowed us to investigate 
the spatial properties of the object. We then performed the spec- 
tral analysis of the regions marked in Fig. 3. To compare both 
data sets, we also analysed the spectra based on the same an- 
nuli as used in the XMM-Newton analysis. This comparison is 
discussed in Sect. 4.1. 



3.3. Spectral models 

We fitted the spectra with the SPEX package (Kaastra et al. 
1996). We employ two different models. The first is the colli- 
sionally ionised equilibrium (CIE) plasma model based on the 
MEKAL code (Mewe et al. 1995), implemented as the cie model 
in SPEX. This model represents an isothermal gas, which is ap- 
propriate for low-statistics data, where we cannot distinguish be- 
tween single and multi-temperature models. For data with higher 
statistics, a differential emission measure (DEM) type model can 
be used (e.g. Kaastra et al. 2004). In this model the emission 
measure Y — J n e nudV (where n e and «e are the electron and 
proton densities and V is the source volume) is a function of tem- 
perature. The wdem model (Kaastra et al. 2004) is a parametri- 
sation of this type, which provides a good empirical description 
of the cores of clusters of galaxies (e.g. de Plaa et al. 2004; 
Kaastra et al. 2004). The model consists of thermal components 
distributed as a truncated power law. The model is characterised 
by 



dY 
df 



cT l/a 





T > T„ 



V T < /3T n 



(1) 



The lower temperature cutoff of the distribution is located at 
/3kT max . This form of the model is chosen in a way that if a 
approaches 0, the model becomes isothermal with temperature 
&7max- For this study, we empirically found that ft = 0.18 yields 
a model which has negligible degeneracy between the maximum 
temperature (kT m - ix ) and the power law slope (I /a). This value 
agrees with earlier work (de Plaa et al. 2006; Werner et al. 2006). 

For a comparison between the outcomes of the (multi- 
temperature) wdem model to single temperature-models, we em- 
ployed the mean temperature, r mean . This value is defined as 



fT§dT 



(2) 



The evaluation of the integral between /3kT max and kT m - dx yields, 
in terms of the model parameters 



kT„ 



(1 + 1/or) (1 -J3 [ ' a+Z ) 
(2+ l/a)(l -/?!/«+!) 



kT„ 



(3) 



We used a neutral absorbing gas model for the determination 
of the interstellar hydrogen column density. Since the hydrogen 
column density value from Hi maps (e.g. Dickey & Lockman 
1990; Kalberla et al. 2005) does not include molecular hydrogen, 
it is possibly lower than the value needed to fit the spectrum. 
Therefore we kept A^h as a free parameter in our fits. 

4. Results 

A summary of the parameters fitted to the data from the 3' 
around the centre of the three clusters (which we will hereafter 
refer to as the "integrated" spectrum) is given in Table 4. The 
nominal statistical uncertainties of the fits are shown, as well as 




Energy (keV) 

Fig. 4. The EPIC-pn spectrum of J1539 (crosses) with subtracted back- 
ground level (histogram) as explained in Sect. 3.1.2 and the fitted model 
(connected line). The part above 2 keV is rebinned with a factor two for 
presentational purposes. 

Table 3. Summary of the fits to the Chandra ACIS spectra of J1539. 



region" 


kT 


Fe abundance 




(keV) 


(proto-solar) 


1 


3.02 ±0.13 


0.9 ± 0.3 


2 


2.7 ± 0.2 


0.41 ± 0.14 


res* 


3.28 ± 0.14 


0.52 ± 0.11 



' The regions are shown in Fig. 3. 

' For this spectrum, regions "1" and "2" are excluded 



(in brackets) the uncertainty associated with the statistical uncer- 
tainty of the subtracted background. In Fig. 5, optical and X-ray 
images of the sources are shown. We also fitted the spectra of 
different annuli employing the same model as for the integrated 
spectrum (Figs. 6 and 7). In the following sections we describe 
the results for each source. 



4.1. J1539 

The integrated spectrum of our brightest source, J 1539, is shown 
in Fig. 4. As shown in the plot, the background is high compared 
to the source (see also Sect. 3.1.2). The effect is noticeable for 
the part of the spectrum above 2.5 keV. We fitted a wdem model 
to the background-subtracted spectrum. The mean temperature 
of the source is 3.03±0.15keV. 

Visual inspection of the Chandra image of J1539 (see Fig. 3) 
shows complexity in the cluster morphology: a compact bubble- 
shaped region and an extended structure in the upper part of it 
stand out from the weaker diffuse emission of the cluster. The 
three regions indicated in the image yield the parameters cited 
in Table 3. Due to the relatively low statistics, we fitted a single- 
temperature model to these spectra. The region labelled "2" has 
a temperature of 2.7+ 0.2 keV, which is slightly below the tem- 
perature of the surrounding gas of 3.28 ± 0.14keV. This is con- 
sistent with the values found in the XMM-Newton analysis. Fits 
to the 0-3' region of the Chandra observation yield a hydrogen 
column density of (17.2 + 1.3)xl0 24 m~ 2 . This is consistent with 
the best-fit column density derived for the XMM-Newton data, 
(20.0 + 0.4 (± 0.2))xl0 24 trT 2 . Here, the uncertainty between 
parentheses is due to the uncertainty in the background deter- 
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Fig. 5. X-ray (top) and optical (bottom) images of the studied clusters. Top and bottom image are at the same scale and the length of the scale 
bars is 30". J1539 and J0812 both have a clear cD galaxy in the centre. In J1601 there are three bright foreground stars in the field of view, which 
makes any optical central structure in this cluster hard to see. 



Table 4. Best-fit parameters for spectra of the 3' around the centre of the galaxy clusters. 



Parameter 


J1539 (XMM-Newton) 


J1539 {Chandra)" 


J1601" 


J0812 


Background scaling factor 


0.99 + 0.02 


0.976 ± 0.013 






Hydrogen column density (10 24 irT 2 ) 


20.0 + 0.4 (± 0.2) 


17.2 ± 1.3 


6.6 ± 0.4 (± 0.3) 


29 ±2 


Hydrogen 21 cm column density^ (10 24 m~ 2 ) 


8.38 


8.38 


6.25 


22.1 


Emission measure^ (10 73 irT 3 ) 


2.57 ± 0.04 (± 0.01) 


2.44 ± 0.07 


9.2 ± 0.2 (± 0.2) 


2.80 ±0.16 


Temperature (keV) c 


3.03+ 0.15 (± 0.08) 


3.50 ±0.11 


8.24 ± 0.58 (± 0.97) 


3.2 + 0.2 


l/af 


0.48 ± 0.09 (± 0.03) 


0.48 (fixed) 






Fe abundance 


0.61 ± 0.03 (± 0.02) 


0.57 ± 0.07 


0.58 + 0.10 (± 0.06) 


0.37 ±0.10 


Si/Fe abundance 


0.65 ± 0.12 (± 0.01) 








Ni/Fe abundance 


2.81 +0.56 (± 0.02) 








Flux (0.1 -2.4 keV) [lO^WnT 2 ] 


1.03 ±0.02 (± 0.01) 


1.03 ±0.03 


0.99 ± 0.02 (± 0.02) 


0.66 + 0.04 


Flux (2-10 keV) [lO" 14 W nr 2 ] 


1.11 ±0.02 (± 0.04) 


1.23 ±0.03 


1.70 + 0.04 (± 0.13) 


0.88 + 0.05 


C-stat/d.o.f. 


270/229 


389/295 


252/224 




* 2 /d.o.f. 








189/134 



The first quoted uncertainty is the statistical uncertainty on the value. The uncertainty between parentheses is due to the uncertainty in 

background estimation (Sect. 3.1.2). 

* For this source, we fitted a "wtfem" model. 

c The emission measure is defined as Y = § n e nu&V, see also Sect. 3.3 
d Deduced from the 21 cm map by Dickey & Lockman (1990). 
e For wdem: "mean temperature", see Eq. (3) 
f As defined in Eq. (1). 



mination by comparing the values of maximum and minimum 4.2. J 1601 
background scaling factors. 

The temperature fitted for the integrated spectrum of this source 
is 8.24+0.58(+0.97)keV (see Table 4). The fit is acceptable (see 
Table 4). We did not find significant multi-temperature structure 
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Fig. 6. Radial profiles of temperature, Fe abundance and Ni/Fe abun- 
dance ratio for J1539. The solid crosses represent the values for sub- 
tracting the background using the vignetting factor listed in Table 4. The 
dashed crosses represent the result when subtracting the background 
using the vignetting factor minus its statistical uncertainty, the dotted 
crosses represent the case when the vignetting plus its statistical uncer- 
tainty is used for subtraction. The dash-dotted diamond shapes represent 
the values fitted to the Chandra ACIS data in the same spatial region. 
The grey stars represent the values fitted by Cavagnolo et al. (2009) us- 
ing the Galactic H i value for the hydrogen column density. The arrows 
in the bottom panel represent 2cr (95.4% confidence level) upper limits. 



for this source. Figure 7 shows that there is a flat radial temper- 
ature and iron abundance distribution. Silicon and nickel abun- 
dances were not measured in any significant way. 



4.3. J0812 

This is the only source for which the background during the ob- 
servation was not high. Because we multiplied the counts in the 
image with the theoretical vignetting weighting, the uncertain- 
ties are not the square root of the measured counts anymore. In 
this case, the Poisson statistic is not valid for the observation, 
and we used a fit based on the x 2 - The temperature of the in- 
tegrated spectral region is measured to be 3.2+0.2keV. The Fe 
abundance is 0.37+0. 10 solar. No other abundances could be sig- 
nificantly determined. 
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Fig. 7. Radial profiles of temperature and Fe abundance for J 1601. The 
solid crosses represent the values for subtracting the background using 
the vignetting factor cited in Table 4. The dashed crosses represent the 
result when subtracting the background using the vignetting factor mi- 
nus its statistical uncertainty, the dotted crosses represent the case when 
the vignetting plus its statistical uncertainty is used for subtraction. 



5. Discussion and Conclusions 

We reported the spectral properties of three moderately bright 
clusters, which were selected from WFC data. Radial tempera- 
ture and abundance profiles of the the cores of J1539 and J1601 
were obtained. In addition, the temperature and iron abundance 
in the highly absorbed cluster J0812 are measured. We were 
able to derive accurate temperatures and abundances, despite the 
high soft-proton background in two of our cluster observations. 
In Sect. 3.1.2, we described a novel method to estimate the vi- 
gnetting factor of soft protons for EPIC-pn. The systematic un- 
certainties introduced by this new method are mostly similar to 
or smaller than the statistical uncertainties found in our spectral 
fits. From the plots in Figs. 6 and 7, it is clear that the effect 
of the background subtraction becomes more important further 
away from the centre of the field of view due to the distribution 
of the source brightness. Especially in the 2-3'region of the up- 
per panel of Fig. 6 it is clear that this effect is non-negligible 
when measuring temperatures. 

The difference we fitted between the 0. 1-2.4 keV fluxes we 
found as well as the fluxes from both the ROSAT all sky sur- 
vey (Voges et al. 1999) and the REFLEX catalogue (Bohringer 
et al. 2004) are typically 20-30%. A comparison of the fluxes 
with catalogue values is however not straightforward. The mea- 
sured flux depends on the model used and the extraction re- 
gion (in our case 3')- Furthermore, the hydrogen column density 
used in these catalogues is based on the H i data from Dickey & 
Lockman (1990). The value we used is found by fitting the spec- 
tra. We found a higher column density in two of our cases. For 
J1539, the higher column density is also found with the Chandra 
data. 
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5.1. J1539 



J1539 is the brightest source of our cluster sample. We found 
a hydrogen column density that is a factor 2.4 higher than the 
value from Dickey & Lockman (1990). This value is consis- 
tent between XMM-Newton and Chandra. The mean temper- 
ature of the integrated spectral region for this cluster is about 
3keV, which indicates that J1539 is a relatively cool cluster. 
This source also had a significant multi-temperature structure 
when we fitted the spectrum with a differential emission mea- 
sure (DEM) model. The strong radial iron gradient in J1539 is 
similar to abundance gradients found in cool-core clusters (De 
Grandi et al. 2004; Tamura et al. 2004). In the Chandra data, no 
clear trend is visible, which is due to the large error bars of the 
outer data points. 

The Chandra image shows structure in the core (see Fig. 3). 
The central elliptic region (labelled "1") is cooler than its sur- 
roundings, and the top region (labelled "2") is cooler than the 
rest of the gas, which is at comparable distances from the core 
(labelled "res."). The Fe abundance is consistent with the value 
derived from the XMM-Newton observation. There is a slight 
hint for a higher Fe abundance in the region "1". The morphol- 
ogy looks similar to other examples of cool-core clusters where 
bubbles were observed, for example in Abell 2052 (Blanton et al. 
2001) and Hydra A (McNamara et al. 2000). The bubbles in 
clusters are thought to be caused by AGN outbursts in the cen- 
tral galaxy (for a review on this topic see McNamara & Nulsen 
2007). These outbursts are detected in the radio domain and cre- 
ate cavities in the X-ray emission of the gas. These AGN also 
seem to provide the energy needed to keep the gas from cool- 
ing below ~ 0.5 keV. Gas at these temperatures is not observed 
in clusters of galaxies (Tamura et al. 2001; Kaastra et al. 2001; 
Peterson et al. 2001, 2003; Kaastra et al. 2004), therefore a re- 
heating process of the gas is needed. 

In the top panel of Fig. 6 we present the temperature profiles 
derived from our fits, and compared them to the profile derived 
by Cavagnolo et al. (2009). Both our radial XMM-Newton and 
Chandra temperature profiles are significantly different from 
their profile because they keep the hydrogen column density 
fixed to the H i value. We let the hydrogen column density free 
in our fits, and found a value which is 2.4 times higher in both 
the Chandra and the XMM-Newton data sets. 

The radial temperature profile derived from the XMM- 
Newton data shows a strong decline in the 1-3' region. The 
Chandra data however show a much weaker decline. We com- 
pared the radial temperature profiles with a profile, based on av- 
erage profiles compiled by Leccardi & Molendi (2008). They 
show profiles for the cooling-core, unidentified cooling core and 
non-cooling core clusters. This classification is based on the 
radial behaviour of temperature in the cluster core. Since we 
did not measure a significant temperature decrease with XMM- 
Newton in this region, we compare our results with the average 
of the unidentified cooling core clusters. 

The temperature kT mean in Leccardi & Molendi (2008) is de- 
fined by fitting a single temperature model to the region between 
0. l/?i8o and 0.6Rm l . However, the annuli we fitted extend only 
to 0.2/?200- Therefore we could not use this range for our mean 
temperature. To compare the two data sets with the trend, we 
scaled the relative temperature to the value of the innermost bin. 
This results in the plot shown in Fig. 8. 



0.05 



0.1 



0.15 



0.2 



R/R 2 , 



1 In Leccardi & Molendi (2008), Risq is used. Based on their defi- 
nition of that radius and the definition of the R 2 oo from Arnaud et al. 
(2005) we find thatfi 20 o = 1.012tfi 80 . 



Fig. 8. Comparison of the XMM-iVewton and Chandra radial tempera- 
ture profiles found with the plot of the uncertain cool-core cluster sam- 
ple from Leccardi & Molendi (2008). The small crosses represent the 
data points from Leccardi & Molendi (2008). The lines represent the 1 
<r upper and lower limit on the sample. The large crosses are the data 
points of our XMM-Newton observation. The crosses with triangles are 
the points obtained by fitting the Chandra data. Because of the differ- 
ence in definition of the mean temperature, the temperature values are 
all scaled to the innermost bin. The R 2 qo is calculated using the relation 
derived by Arnaud et al. (2005) 



Since the shape of the Chandra and XMM-Newton temper- 
ature profiles are different, the scaling on the innermost bin re- 
sults in more significant deviations between the profiles at larger 
radii (see Fig. 8). The absolute temperatures in the three inner 
bins lie within 2.1<x. Furthermore, the scaled XMM-Newton data 
point at 0.17/?200 is 4.7cr below the average profile of Leccardi 
& Molendi (2008). 

The difference in profile shapes between XMM-Newton, 
Chandra and the average profile can be explained in multiple 
ways. The cooler temperature of the innermost bin as measured 
by Chandra may be partly due to point spread function (PSF) ef- 
fects. The larger PSF of XMM-Newton may cause photons from 
the second radial bin to scatter into the central bin and subse- 
quently raise the derived temperature. In addition, the subtracted 
soft-proton background may not exactly represent the true soft- 
proton spectrum as measured in the annuli. This is because the 
soft-proton spectrum in the inner annuli might have a slightly 
different slope compared to the soft-proton spectrum in the out- 
ermost annulus on which the subtracted background is based. 
The slope of the subtracted background can have a significant 
influence on the measured temperature, especially in the outer 
annuli where the background is high with respect to the source 
spectra. 

5.2. J1601 

The hydrogen column density that we fitted for J 1601 is consis- 
tent with the value based on the H i maps of Dickey & Lockman 
(1990). We derived a temperature for the integrated spectrum, 
which is rather high for clusters of galaxies. Clusters with high 
temperatures have less prominent spectral lines, and the fit of the 
temperature is mostly dependent on the continuum. Therefore, 
the uncertainty on the background scaling factor introduces a 
large systematic uncertainty on the temperature. 
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The flat temperature and iron abundance profiles suggest 
however that J 1601 is a non-cooling core cluster. This means ei- 
ther that the core of the cluster has not cooled yet or was recently 
reheated. The rate at which the cluster core cools can be esti- 
mated by the cooling time. This is the typical time scale during 
which the core should have radiated all its energy away. Sarazin 
(1988) gives for the cooling time of a cluster of galaxies 

^= 8 - 5xi °iTo^r(io^) 1/2yr - w 

In this equation, T g is the gas temperature and n p represents the 
proton density. The volume is calculated with the Cosmology 
Calculator by Wright (2006). We used the emission mea- 
sure that we found for the innermost 30" of the cluster core 
((1.98 + 0.08) x 10 73 irT 3 ) and the temperature of the same re- 
gion (7.4+0.7 keV), assuming that this region is a sphere with 
constant density. We took the electron density to be a factor 1 .2 
higher than the proton density. The cooling time we compute is 
4.9+0.2 Gyr, which is about 1/2 of the Hubble time at the red- 
shift of J 1601, which indicates that the cluster has had enough 
time to cool. 

A merger at an early stage of the clusters life could cause 
shocks, which would reheat the gas. This could explain the ab- 
sence of a cool core when the cooling time is more than ~ 0.1 
Hubble times (e.g. Burns et al. 1997; Gomez et al. 2002). The 
image of this cluster made using XMM-Newton (top-middle 
panel of Fig. 5) and the intensity contours of the ROSAT HRI 
(Fig. 1) hint towards an elongation of the X-ray cluster emission. 
This hint is however not statistically strong enough to pinpoint 
a merger. The bottom-middle panel of Fig. 5 shows three very 
bright foreground stars, which make it hard to state whether there 
is one central galaxy or if the core consists of multiple galaxies, 
which could indicate a merger history. 

5.3. J0812 

The 0. 1-2.4 keV flux we derived for this source is 60% lower 
than the ROSAT value of Voges et al. (1999). The integrated 
temperature is relatively low and the iron abundance is slightly 
lower than the average of larger cluster samples (e.g. Tamura 
et al. 2004; De Grandi et al. 2004; de Plaa et al. 2006). The 
hydrogen column density we fitted is high, 36% higher than the 
value of Dickey & Lockman (1990). Although this is the only 
observation for which the soft-proton background is very low, it 
is also the dimmest of the three sources. Therefore, the statistical 
quality of this observation is limited. 
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